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lated with the three-loop accuracy. The last result can be used to improve the accuracy 
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coefficient for the flavour-nonsinglet QCD current with n antisymmetrized 7-matrices is 
also obtained at three loops; the result for the tensor current (n = 2) is new. 
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1. Introduction 

Let us consider QCD with a heavy flavour, say c. It is well-known that processes with light 
quarks and gluons having characteristic momenta much less than m c can be described by 
an effective low-energy theory — QCD without c-quarks. The renormalized light-quark 
and gluon fields A(fj,) in the full theory are related to the corresponding fields in 

the effective theory, up to corrections suppressed by powers of l/m c , by the decoupling 
relations 

«(//) = C} /2 (^M , A<ji) = cYWM , (1.1) 

where all quantities in the low-energy theory are denoted by primes (we use the MS renor- 
malization scheme throughout this paper). Similarly, the coupling constant, the gauge- 
fixing parameter and the light-quark masses in the two theories are related by 

"sO) = Ca(A*)«a(A*) > a (^) = (,A{p)a'{n) , m i{^) = CmivWiiv) • (1-2) 

The QCD decoupling coefficients are known at three Q and even four loops [§, § 1 . Various 
operators of full QCD, e.g., light-light quark currents, can be expressed via operators of 
the low-energy effective theory, similarly to (|1 . 1| ) . 

Now let us consider the 6-quark HQET. If the characteristic residual momentum of 
b, as well as characteristic momenta of light quarks and gluons, are much less than m c , 
then all processes can be described by the low-energy effective theory — HQET without 
c-quarks. The decoupling relations for the light fields ( |1.1|) and the parameters of the 
Lagrangian ( |1.2| ) are exactly the same as in QCD. The static 6-quark field Q in the "full" 
theory (HQET with c-quarks) and in the effective theory (HQET without c-quarks) are 
related by 

Q(m) = cJ /2 (m)Q'(m)- (1-3) 

Various operators of the "full" HQET can be expressed via operators of the low-energy 
HQET. For example, for the heavy-light quark currents j = qTQ we have 

j» = C>)iV), (1.4) 

up to corrections suppressed by powers of l/m c . The HQET-specific decoupling coefficients 
Cq, Cj have been calculated in || with the two- loop accuracy 2 . Decoupling in HQET is 
also discussed in ||, Sects. 4.7 and 5.5. 

In the present paper, we shall calculate Cq and Q up to three loops. To do so, we need to 
calculate on-shell HQET propagator diagrams containing a massive quark loop. Reduction 
of scalar Feynman integrals of this type to master integrals is considered in Sect. [| the 
master integrals are calculated in Sect. ||. The HQET decoupling coefficients are obtained 
in Sects. ^| and ||. In Appendix |A| we derive the three-loop decoupling coefficient for the 
QCD flavour-nonsinglet quark current with n antisymmetrized 7-matrices, as a generic 
formula containing n; the result for the tensor current (n = 2) is new. 

lr The result of || contains one master integral which was not known analytically, only numerically, with 
37-digits precision. An analytical expression for this integral has been published later [^] . 

2 One more HQET decoupling coefficient, that for the 6-quark chromomagnetic interaction, has been 
found in M, based on the calculations of M, also at two loops. 
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2. Reduction of Feynman integrals 



2.1 General remarks 



In this subsection, we discuss some general prop- 
erties and simple particular cases of on-shell HQET 
integrals with a massive quark loop. 

Let us consider the integrals (Fig. [l]) 



F (0,1,0,2) = J -^aT 
Ex = -2k -v-iO, 



f{k 2 )d d k 




rpa 2 ) 
^2 

En 



Figure 1: Diagram with a single 



iO (2.1) HQET line 



with some f(k 2 ). Following the strategy of integration 
by parts || and integrating the identity 



d_ 

dk 



E x ) E^E^ 2 
we obtain the recurrence relation 



d-at-2- 4(ai + 1) 



E2 

eJ 



f(k 2 



E^E^ 2 



(2.2) 



(d - ai - 2)F(ai,a 2 ) = 4(ai + 2~ F( ai ,a 2 ) 



(2.3) 



(here, as usual, the operator 2 decreases a 2 by 1, and 1 ++ increases a\ by 2). Its solution 
is 



F(a 1 ,a 2 ) 



-4) 



-ai/2 



r(f) iW) 



' d— ai 



«1 



r(^±i)r(^) 



> -F(0,a 2 + yj , even 01, 

ai-l\ (2.4) 
F(l,a 2 + -^— J , oddai >0, y ' 

odd ai < 0, 



The result for even ai < can easily be derived by averaging over v directions in the 
d-dimensional Euclidean space after the Wick rotation; it was supposed in |l0| that the 
same formula can be used for even a\ > 0, and this method was used, e.g., in 0, ||, but 
no proof existed. Now we see that this averaging trick follows from ( |2.3| ). For odd a\ < 0, 
averaging over v directions gives 0, but this result does not extend to odd a\ > 0. 
Now we shall discuss the two-loop integrals (Fig. §) 



F(ai,a 2 , 03,04,05) 



d d h d d k 2 



h + k 2 



(zvr rf / 2 ) 2 J El x EfEfE a *Ef 

(2.5) 



with four denominators 




Ei 
E 3 



—2k\ ■ v — iO . 
1 - fc| - iO , 



E<> 



-k\ - iO , 



E 4 = l-(ki + h) 2 



HO 



(2.6) 



Figure 2: Two- loop diagram 
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(we have put the quark mass m = 1; the power of m 

can easily be restored by dimensionality), and one numerator 



E 5 = (2k 2 + ki) ■ v (2.7) 
(05 is always < 0). Integrals with 03 < or 04 < vanish. We have 

F(ai, 0,2,0,3, 0,4,0,5) = (-l) a5 F(oi, 02,04,03,05) . (2.8) 

These integrals can be subdivided into two disjoint subsets, with even and odd a\ + 05. 
We shall call them "apparently even" and "apparently odd" (they would be even and odd 
with respect to v — * —v if there were no — iO in E\). Apparently odd integrals with a\ < 
vanish, because we can omit — iO in the numerator. These two classes are not mixed by 
any recurrence relations, therefore, we have two disjoint problems. 

These integrals can be calculated as follows. First, we get rid of the numerator E^ 0,5 . 
We have a contraction of —05 vectors v with the one-loop tensor integral which depends on 
k\. Decomposing this integral into all possible tensor structures, we obtain an expression 
containing only E\ (no E§). Using a-parametrization, we obtain 

F(oi, 03,03, 04,0) = (2.9) 



r (a) r [Si^i + a2 + o 3 ) r + a 2 + 04) r(f + a 2 + a 3 + a 4 - d)T ( ^ 



0,2 



2r(ai)r(a 3 )r(a 4 )r(ai + 2a 2 + 03 + 04 - d)T 



d—a. 



In full accordance with ( |2.4| ), integrals F(a\, 02, 03, 04, 0) with even a\ reduce to F(0,a2 + 
01/2,03,04,0) (this is a well-known two-loop vacuum integral fll|]); those with odd ai > 
reduce to F(l, 02 + (ai — l)/2, 03, 04, 0); and those with odd a\ < vanish. This is a strong 
check. All apparently even integrals are proportional to the single master integral 



^0 = 4^, (2-10) 




where 

is the one-loop vacuum integral. All apparently odd integrals are proportional to the single 
master integral 

jo = -^^V = 2«-v y-M) . (2 . 12) 

u \~~/ r 2 (2-d/2) v ' 

In order to perform our decoupling calculation, we need two classes of three-loop in- 
tegrals. They will be reduced to master integrals using integration by parts S. A similar 



problem (off-shell HQET integrals without massive quarks) has been solved earlier [12] 
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However, constructing reduction algorithms by hand is tedious and time-consuming. Sev- 
eral algorithmic methods to solve reduction problems for Feynman integrals have been 



recently suggested; they are discussed, e.g., in [13]. Some approaches are based on the use 
of Grobner bases Jl4| , |l5| [Iq ]. This is the last algorithm in this list [15, 16] that we have 
used in the present calculations. 

2.2 Integrals of class 1 




ki k 2 
Figure 3: Diagram of class 1 

These integrals (Fig. ||), 

Fi(ai, a 2 , as, 04, a^,aQ,a7, as,ag 



1 



( i7r d/2)E 



d d ki d d k 2 d d k 3 



(2.13) 



have 8 denominators, 



-2ki ■ v , 



Eo. 



-2k 2 ■ v , 



-k 



1 ! 



Ea 



-k 2 



e 5 = -(h - k 2 y 

\2 



(2.14) 



Eh. 
Ez 

E 6 = l-k%, E 7 = l-(kx + k 3 ) 2 , E 8 = 1 - (k 2 + k 3 ) 2 

(here and in what follows, — iO is implied in all denominators), and one numerator, 

E 9 = 2k 3 -v (2.15) 

(ag is always < 0). The integral ( [2.13D is symmetric with respect to (1 *-* 2, 3 «-> 4, 7 «-> 8). 
It vanishes if indices of the following subsets of lines are non-positive: {5, 7}, {5, 8}, {6, 7}, 
{6,8}, {7,8}, {3,4,6}. 



We applied the method of [15, 16] to construct an algorithm for reducing integrals of 
this class to master integrals. Recurrence relations do not mix apparently even integrals 
(with even ax + a 2 + ag) with apparently odd ones. All apparently even integrals reduce 
to linear combinations of the following master integrals: 




h 




h 




(2.16) 



All apparently odd integrals reduce to the following master ones: 



-Wo 




2.3 Integrals of class 2 



h + h 




k 2 + k 3 



Figure 4: Diagram of class 2 



These integrals (Fig. f|), 

F 2 (ai, a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 , a 9 ) 



(ivrrf/2)S 



<rfci d d fc 2 d% 



have 7 denominators, 



£•0 



3 ! 



E* = l-k 



1 J 



a = 1 - fe 



2 > 



E 5 = 1 - (h + k 3 ) 2 , £ 6 = 1 - (fc 2 + fc 3 ) 2 , E 7 = -(k 1 -k 2 ) 2 
and two numerators, 

E 8 = (2ki +k 3 )-v, E 9 = (2k 2 + k 3 )-v 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(ag and a 9 are always < 0). The integral ( 2.18j ) is symmetric with respect to (3 <-> 4, 5 «-> 
6, 8 «-> 9), and 

F 2 (oi, 02,03, 04,05,06, 07,08,09) = (-l) as+a9 F 2 (ai, o 2 , 05,06,03,04,07,08,09) . (2.21) 

It vanishes if indices of the following subsets of lines are non-positive: {3, 5}, {4, 6}, {3, 4}, 
{5,6}. 



We have used the method of ]15| , 16 1 to construct an algorithm reducing these integrals 
to master ones. For apparently even integrals (with even ai + as + ag) we also used a 
method similar to subsection |2.l| . Namely, we have a contraction of —as — a 9 vectors 
v with a two-loop tensor integral which depends on k 3 . Decomposing this integral into 
all possible tensor structures, we get rid of the numerator E^ aa Eg ag . Then we use the 



averaging formula (2.4), and the problem reduces to three-loop massive vacuum integrals 
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B M [y]. A REDUCE package RECURSOR |j reduces them to two master integrals, i, 



(see (|2l6D ) and 

(2.22) 




Agreement of results produced by these two ways serves as a strong check. 

Apparently odd integrals of this class are expressed via IqJq, J\ (see ( 2.17] )), and 



J 3 = /A \ . (2.23) 



3. Master integrals 

3.1 Simple master integrals 

Simplest master integrals can exactly be expressed via V functions: 
T{d/2 - l)r(l - d/2)r 2 (3 - d)T{4 - 3d/2) 



r(6 - 2d) 

m-17 2 r(d/2-i)r(5-2d)r(7-3d) 



J\ = Z IT 



T(2 - d/2)r(3 - d)r(4 - 3d/2) 



7rV^(7/2 - d)T(2 - d/2)T\d/2 - l)r((3(3 - d))/2)r(d - 5/2) 

2 2r(7-2d)r(d-2)r((d-i)/2) ' 

The vacuum integral ^ ( |2 .22[ ) has been investigated in |17|| . It can be written as 
T _ A (d-2) 2 (3d-ll) r g ? d-4 fj 



(d-3)(3d-8)(3d- 10) u 2d - 7 

+ 64 ( rf -2)(d- < 3 K3d 4 -8)(3d-10) r3 " + ^ (£) ' <3 ' 2) 
where d = 4 — 2e, B±(e) can be expressed via 3-F2 of unit argument, and 

B A {e) = B 4 + 0(e), 

B 4 = 16 Li 4 Q) + I In 2 2 (In 2 2 - . 2 ) - J|vr 4 ; (3.3) 

its expansion up to the e 3 term was presented in JT^| . 
3.2 Integral I 2 

We evaluated all the complicated master integrals (i.e. which are not expressed in terms 
of gamma functions for general d) by the method of MB representation which is based on 
the following formula 

1 11 r+ioo -yz 

' d zT(X + z)T(-z)— -. (3.4) 



(X + YY T(A) 2vrz 7„ JOO v ; v > X*+* 
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Here the contour of integration is chosen in the standard way: the poles with a T(. . .+z) 
dependence (let us call them left poles, for brevity) are to the left of the contour and the 
poles with a T(. . . — z) dependence {right poles) are to the right of it. This formula is used 
to replace a sum of terms raised to some power by their products raised to some powers at 
the cost of having an extra integration. 

Two systematic strategies for evaluating MB integrals in a Laurent expansion in e were 
suggested in [19, 20]. A description of this method is presented in [21] and Chap. 4 of [13]. 

The simplest variant of using (|3.4] ) is to write down a massive propagator in terms of 
massless ones [22|: 

1 11 r+ io ° (m 2 Y 

1 d. \> r(\ + z)r(-z). (3.5) 



(m 2 - k 2 ) x T(A) 2m J_ ioo (-k 2 ) x + z 

To have more checks, it is very useful to derive such MB representations for general 
indices. To evaluate I 2 we derived a MB representation for Fi(oi, 02, 0, 0, 0, Og, 07, ag, 0). 
We applied ( |3.5| ) to the last two factors in the integrand (with 07 and a%) and evaluated 
the two resulting integrals over k\ and k 2 and the resulting integral over kj, by well-known 
one loop integration formulae to obtain 

Fi(a 1 ,a 2 ,0, 0, 0, a 6 , a 7 , a 8 , 0) 



x (0 -\2 / / dzidz 2 r(oi 2 /2 + a 678 + 3s - 6 + z\ + z 2 ) 

[ZTTl) J —too J —ioo 

r(ai/2 + a 7 + £ - 2 + zi)T(a 2 /2 + a 8 + e - 2 + z 2 ) 
r((oi + a 2 - 7)/2 + a 7 + a 8 + 2e + z 1 + z 2 ) 
xT((ai - 3)/2 + a 7 + s + zi)T{{a 2 - 3)/2 + a 8 + e + z 2 ) 

xr(2 - a 7 - £ - zi)T(2 -a 8 -e- z 2 )T(-z 1 )T{-z 2 ) , (3.6) 

where a\ 2 = a\ + a 2 etc. 

It was more convenient to obtain a Laurent expansion of 1,0, 0, 0, 1, 1, 1,0) in 
e from a Laurent expansion of the integral ^(2,2,0,0,0,1,2,2,0) and its reduction, in 
particular, to I 2 . For this auxiliary integral, we used (|3.6|) . The resulting twofold MB 
integral can be expanded immediately in e because there is no gluing of poles of different 
nature when e — ► 0. Then one can close the integration contours over z\ and z 2 to the 
right and obtain a double series. Its summation gives the following result: 

Fx (2, 2, 0, 0, 0, 1, 2, 2, 0) = T 3 (l + e) - ± - + 0{e 2 )^j (3.7) 

which leads to 

/ 2 = -r 3 (l + e )^(i + ^ + 0(e)) . (3.8) 

Integrals Fi(ai, a 2 , 0, 0, 0, ae, a 7 , ag, 0) can be calculated also in the coordinate space. 
After continuation to Euclidean time, we have 

2<5d/2— a\ —a2—aQ—a7—as+3 

Fi(ai, O2,0, 0,0, 016,07,08,0) 



r(a 1 )r(a 2 )r(a 6 )r(a 7 )r( c 
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f'OO poo 

Jo Jo 

K d/2 -a 7 {ti)K d/2 ^ as {t 2 )K d/2 _ a6 (h + t 2 ) . (3.9) 

If we expand the Bessel functions K d / 2 _ a7 (ti) and K d / 2 _ a8 (t 2 ) in ti and t 2 , then the 
integrals can be calculated, and we obtain a double series. However, this series is not very 
convenient for expansion in s, because even for convergent integrals separate terms contain 
1/e 2 . We have verified that the 0(1) term in ( p. 7] ) is reproduced by numerical integration 
in (|3.9| ) with 12 digits accuracy, thus providing a good check. 



3.3 Integral I 3 

To evaluate I3 we derived a general MB representation for F\(a\, a 2 , 03, 04, 0, 0, 07, ag, 0), 
where ^3 = i*i(l, 1, 1, 1, 0, 0, 1, 1,0). To do this, we used an alpha representation and then 
applied (|3.4p twice in an appropriate way. We obtained 



„, n n n x r(oi 2 /2 + a 347 + 2e - 4)r(a 12 /2 + a 348 + 2e - 4) 

Fi (01,02,03,04,0,0,07,08,0) = ^-p//, q \ T^yrV To 0\ 

4vvrr((4 - 012 - 2e)/2)r(oi2 + 2a 3 47 8 + 4e - 8) 
£(012/2 + Q3478 + 3e - 6) 1 /- +lo ° r(a 1 /2 + z 2 )r(l/2 + a 1 /2 + z 2 ) 

nr(a.) (2^) 2 Loo 7_ ioo Zl 22 r(i/2-*!) 

xT(ai 2 /2 + 2fi + z 2 )r(a 34 + e - 2 - zj)r(2 - ai 2 /2 - a 4 - e - z 2 ) 

r(-oi/2 - zi - z 2 )r(l/2 - oi/2 - zi - z 2 )r(2 - a 3 - e + zi + z 2 )F(-z 2 ) . (3.10) 

As in the previous case, we preferred to calculate, instead of I3, another integral of 
this family, i*i(2, 2, 1, 1, 0, 0, 1, 1,0) and then obtain ^3 using our reduction procedure. We 
have, after some changes of variables, 

r(3e)r(l + 2e) 2 

Fi(2, 2, 1,1, 0,0, 1,1,0) - 



40Fr(-e)r(2 + 4e) 

+ ' x , . r(3/2 + z!)r(i + z 2 )r(3/2 + z 2 )r(i + zi) 

ClZlClZ 2 



(27Ti) 2 7_ loo J_ loo r(5/2 + z 1+ z 2 ) 

xr(2 + e + zi + z 2 )T(-l -e- zi)T(-zi)r(-l - e - z 2 )T(-z 2 ) . (3.11) 

After the replacement 

r(2 + e + z 1 + z 2 ) -> (1 + e + Zl + z 2 )r(l + e + z x + z 2 ) 

the integral can be decomposed into two integrals where the integration can be performed 
using the first Barnes lemma. Resulting onefold integrals can be also evaluated with the 
help of the first and the second Barnes lemmas so that we obtain a result in terms of 
gamma functions for general d: 

= r(i/2 - e)r(-g)r(2g) 2 r(i + e)r(3e - 1) 

3 ~ 4r(3/2 - e)r(4e) 

x [V>(l/2 - e) + V(l - e) + 2 In 2 + 2 7E ] . (3.12) 
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3.4 Integral I4 

To obtain a MB representation for the integrals Fi(di, 02, 0, 04, 05, og, 07, 0, 0), with I4 = 
.Fi(l, 1, 0, 1, 1, 1, 1, 0, 0), we replace the internal integral over k 2 by a onefold MB integral, 
using Feynman parameters straightforwardly. After this the resulting integral over k\ and 
&3 can be evaluated in terms of gamma functions and we obtain a general onefold MB 
representation which gives 

_ r(2 £ ) 2 r(3e - 1) 1 f +io ° T(l + z)T(l/2 + e + z)T(l + e + z) 



4T(4e) 2m J _ ioo T(3/2 + e + z)T{l - 2e - z) 

xT(-2e - z)T(-e - z)T(-z) . (3.13) 

This integral can straightforwardly be evaluated by expanding the integrand in e, closing 
the integration contour to the right and summing up resulting series. We obtain 



h = -r 3 (i + e 



f_ _ 6C3 - 5tt 2 U_ 4 _ 10 19 2 
9e 2 9e + 270^ 3 ^ 3 + 9 * 

(3.14) 



i 8 . 8 9 . ll/i 38 „ 65 r, \ .) 
+ |-3 &+ 9' Cs+ 54' 4 -S 6 + T r ) 6 + 0(6 



3.5 Integral J3 



To obtain a MB representation for the integrals F 2 (ai, a 2 , 03, 04, 05, a^, 0, 0, 0), with J3 = 
^2(1,0,1,1,1,1,0,0,0), we use the following onefold MB representation for a self-energy 
one-loop integral with two equal masses which can straightforwardly be derived using 
Feynman parameters: 

d d k 1 



m d l 2 J {-k 2 + m 2 )^ [-(q - k) 2 + m 2 } a ^ T(ai)T(a 2 )(m 2 )^+ a ^- 2 

2TiJ- ioo U ! j r(o 1 +(i 2 + 22) 

Writing down the subintegral over k 2 using ( 3.15D we obtain an integral which can 
explicitly be evaluated in terms of gamma functions for general e so that we obtain only 
onefold MB representation for the integral under consideration: 

r(a x /2) 



F 2 (ai, 02, 03, 04, 05, a%, 0, 0, 0) 



2\[ l _ t2 T{ai)T{2 - ai /2 -e)(m 2 ) a i/ 2+a 23456-6+3 £ 



,+ 4 oo r(fl4 + z)r(ag + z )r(_ z )r( 0l /2 + a 23 5 + 2e - 4 - z) 



2ni J_ loo ' r(a 46 + 2z)r(ai 2235 + 2e - 4 - 2z) 

xT(2 - ai/2 -a 2 -e + z)T(a m + e - 2 + z) 

xr(oi/2 + a 23 + e - 2 - z)T{ ai /2 + a 25 + e - 2 - z) . (3.16) 



This gives 



n 3/2 1 r+ioo T{ i + z ^ v ^i2- £ + z )Y(e + z) 
J3 = ■ — ,~ r- — : / dz 



4 e r(3/2 - e) 2m J_ ioo T(3/2 + z)T(e - z) 

xT(-l/2 + £- z)r(-3/2 + 2e - z)T(-z) . (3.17) 
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For our purposes, it was sufficient to evaluate only the leading term of this integral (using 
the strategy of [|l9|| ) in expansion in e: 



397r 2 

J 3 = -^- + 0(e). (3.18) 

4. Decoupling for the heavy-quark field 
4.1 General formulae 

It is convenient to find the relation between the bare fields in the "full " and the low-energy 
HQET, 

Qo=(C°qY /2 Qo, (4.1) 



first; then 



Z' n (a' (a), a' (a)) 
( Q (n) = $± S ^H • (4.2) 

Z Q (a s {ii),a{n)) 



The renormalization constant Zq can be reconstructed from the three-loop anomalous 
dimension of the HQET heavy-quark field Q which has been found in [^jj by an on-shell 
massive QCD calculation and confirmed by an independent HQET calculation in p4|) . 

The bare heavy-quark propagators in the two theories near the mass shell (u — > 0) 
behave as 

S(u) = -2- + . . . , S'(lo) = -9- + ..., (4.3) 

where Zq and Z'q are the on-shell renormalization constants of the heavy-quark field in 
these theories. Therefore, 

Z os 

<7/os ' (4-4) 
The bare heavy-quark propagator can be written as 

S(u) = 1— , (4.5) 

OJ — 1j{UJ) 

where — is the sum of all one-particle- irreducible self-energy diagrams. Therefore, 

Z° Q S = , , (4.6) 

1 - (dE(u)/doj) 



uj=0 



and similarly 



1 - (dfl'(u)/dw) 



(4.7) 



oj=0 

All diagrams for (dE' '(u) / 1 600)^=0 contain no scale and hence vanish, and 



Z'8 S = 1. (4.8) 
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4.2 Bare calculation in full HQET 

Only diagrams with (at least one) c-quark loop contribute to (dEi(u) / (kS) w =Q. Such dia- 
grams first appear at two loops (Fig. ||a). This diagram has been calculated in [10|: 

/ ~ \ (2) 
/ dS(, 



\ dw 



g^m-t r2 (d-l)(d-2)(d-6) 
jr — — — r^-in 



uj=0 



(4vr) c 



'0' 



2(d-5)(d-7) 



(4.9) 



(of course, go and the bare c-quark mass m c o are those of the full theory) . 






+s 




+s 







m 



Figure 5: Diagrams for (dil{uj)/duj)u=Q 

Three-loop diagrams contributing to (cE(uj) / duj)u=o are shown in Fig. ^b-m, where the 
shaded blob means the sum of the massless one-loop insertions (light quark, gluon, ghost), 
and "+s" means that there is also a mirror-symmetric diagram. The diagram Fig. ||i has 
three HQET denominators depending on only two loop momenta; these denominators are 
linearly dependent, and we can eliminate one HQET line by taking partial fractions. After 
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that, all diagrams of Fig. [B|b-j can be expressed via scalar integrals of class 1 (Sect. 
those of Fig. ||k-m can be expressed via scalar integrals of class 2 (Sect. ^ ). 
The three-loop result can be written in the form 

l "ST J - CFTF UM2^2^ CiI i s «> (410) 

\ / lv=0 v ' i=l j=l 

with the colour structures 

Ci=C7 F , C 2 = C A , C 3 = T Fni , C A = T F , (4.11) 



and the integral structures 



{d -2) 



j _ v" -i t3 

1 24(d - 3)(d - 4) 2 (d - 5) 3 (d - 6)(d - 7)(d - 8)(d - 10) ' 

- _ (3d-8)(3d- 10) 

2 ~ 16(d - 3)(d - 4)(d - 5) 2 (2d - 7)(2d - 9)(2d - 11) U 

- 3d - 10 - 3d - 10 

J3 = tt — — : — — : rrris , i/t 



(d-5)(2d-9)(2d- 11) °' (d-5)(2d-9)(2d- 11) ' 

- _ (3d -8) (3d - 10) 

5 ~ 16(d-4) 2 (d-6)(d-8)(d- 10) 5 ' ( ' 12) 

All the non-zero coefficients are: 

s n = 6(d - 5) 2 (d - 10) (2d 8 - 75d 7 + 1212d 6 - 11042d 5 + 62070d 4 - 220131d 3 

+ 478504d 2 - 576108d + 287408) , 
812 = 2(d - 2)(d - 3)(d - 5) 2 (2d - ll)(2d 4 - 29d 3 + 157d 2 - 366d + 288) , 
s 15 = -2s 25 = -2(d - 10) (d 4 - 20d 3 + 141d 2 - 398d + 328) , 
s 21 = (d - 7)(d - 10)(15d 8 - 555d 7 + 8858d 6 - 79542d 5 + 438503d 4 - 1514103d 3 

+ 3179372d 2 - 3678156d + 1767408) , 
S22 = -(d - 3) (2d 8 - 63d 7 + 864d 6 - 6725d 5 + 32377d 4 - 98156d 3 + 181403d 2 

- 184490d + 77536) + 2(d - l)(d - 3) 2 (d - 4)(d - 5) 3 (1 - o ) , 
s 23 = -S24 = 2(2d 2 - 19d + 47) , s 32 = 16(d - 2)(d - 3)(d - 4)(d - 5) 3 , 
s 4 i = 12(d - 4)(d - 5) 2 (9d 6 - 277d 5 + 3387d 4 - 20943d 3 + 68428d 2 - 110236d + 66352) , 
S45 = 4(d - 4)(d 3 - 15d 2 + 54d - 8) . (4.13) 

The coefficient of Cjyli is not gauge invariant. Therefore, the on-shell renormalization 
constant of the HQET heavy-quark field (^~J) is not gauge-invariant, starting from three 



loops. The same phenomenon has been observed in QCD 23]. In the abelian case, the 
on-shell renormalization constant is gauge invariant to all orders of perturbation theory; 
the gauge dependence only appears in the non-abelian colour structure Ca- 



- 13 - 



4.3 Renormalized decoupling coefficient 



Now we calculate the renormalized decoupling coefficient ( |4.2|) . To this end, we re-express 
Zq via <7q) oo; Z'q via g' 2 , a' , which are re-expressed via ao using the bare decoupling 
relations; combine all factors in (p~I^); and, finally, re-express the result via a s (fi), a(fi) and 
expand the coefficients in e. We arrive at 



CqOO = 1 + CpTp I 2L — —L + — 



16 



52 



+ C F T F < 



9 
16 



47T 



T F (nj + 2) 



+ 



+ 



-18C F - i - 3l) C A + 

/121 155\ 16 

6(8C 3 - l)Cjr + ( —a(ji) - 48C 3 + — J C A - -T F (36m + 31) 



+ ( 165 4 - ^ 4 + IO6C3 - C F 



+ 



1 / . 2387\ , . D _ 52 4 413 . 262 

- 8C3 a(u) - 8S4 H vr 4 C3 H 

3 V 72 y VFV 45 9 243 



32 / 179 \ m 8 f 395 \ , 



4tt 



+ 



L 



(4.14) 



where 



L = 2 In ■ 



m c (/i) 

(m c (/Lt) is the MS renormalized mass), and -B4 is given by (|3.3|). 



(4.15) 



5. Decoupling for the heavy— light current 
5.1 General formulae 

The bare heavy-light current Jo = Qo^Qo in the "full" HQET is related to the corresponding 
operator j' in the low-energy theory by 



Jo = 0o 



up to l/m c corrections. The renormalized operators are related by (1.4), where 



(5.1) 



(5.2) 



The renormalization constant Zj can be reconstructed from the three-loop anomalous di- 
mension of the HQET current J which has been calculated in p4| . 

The Green function of jo, Q , and qo can be written as S(p)T(p, u>)TS(u) (Fig. |TJ), 
where T(p,u)T is the proper vertex — the sum of all one-particle-irreducible diagrams 
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(here T is the Dirac matrix in the current; due to the HQET Feynman rules, no 7-matrices 
can appear to the right from T). The Green functions in the two theories are related by 



S(p)f(p,u;)r5H 
Recalling the relations 



1/2 „ 



p,u 



S(P) = (° q S'(p) 



S(u>) = (° Q S>(u) 



(up to l/m c corrections), we obtain 
f(p,uj) 



TTl, 



It is most convenient to use the point p = 0, lo = 0, then power 
corrections may be omitted. The vertices T(0, 0), r'(0, 0) can have 
Dirac structures 1 and f; if the light quark q is massless, the num- 
ber of 7-matrices in them is always even, and the second structure 
does not appear. So, these vertices are scalar, and 



(5.3) 



(5.4) 



(5.5) 




Figure 6: Green func- 
tion 



y/2 f(o,o) 



r'(o,o) 

All loop corrections to r'(0, 0) contain no scale and hence vanish: 

f (0,0) = 1. 

Therefore, 



z5 s z° s V /2 f(o,o). 



(5.6) 



(5.7) 



(5.8) 



^ 1 J 

The right-hand side is an on-shell matrix element of the gauge- invariant operator Jo> an d 
hence is gauge-invariant; the quantities T(0, 0), Zq, Z° s taken separately are not gauge- 
invariant starting from three loops. 

5.2 Bare calculation in full HQET 

Only diagrams with (at least one) c-quark loop con- 
tribute to A(0, 0) = r(0, 0) — 1. Such diagrams first appear 
at two loops (Fig. |7|). As we know, the result is scalar; so, 
we can take j Tr of the 7-matrices associated with the light- 
quark line. In this diagram, there is only one light-quark 
propagator S(k) and the vertex, and taking the trace gives 
k 1 *. The c-quark loop inserted into the gluon propagator is 
transverse, and the diagram vanishes |i~0|| : 




Figure 7: Diagram for A(0, 0) 



A(2)(0,0) = 0. 



(5.9) 
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Three-loop diagrams for A(0, 0) can be obtained from those in Fig. ||b-m by inserting 
the heavy-to-light vertex into all possible places along the heavy-quark line; of course, 
mirror-symmetric diagrams denoted by "+s" in Fig. |5] should be considered separately. A 
number of diagrams with a single light-quark propagator vanish for the same reason as 
the two-loop diagram (Fig. 0). For some diagrams with three HQET lines, we use partial 
fractioning. After that, all diagrams produced from those of Fig. ||b-j are expressed via 
scalar integrals of class 1 (Sect. |2.2| ); diagrams produced from Fig. [|k-m are expressed via 
scalar integrals of class 2 (Sect. |2.3|) . 

The three-loop result can be written in the form 



6 -6s 2 4 
\ ' i=l j=2 



y~i CiijKj 



(5.10) 



(see (lilll) and (jig)). A11 the 

non-zero coefficients are: 



A 12 = -2(d - 2) 2 (d - 4) 2 (d - 5) 3 , Ai 3 = Au = 2(d - 2){d - 5) 2 , 

A 22 = {d- 3) 2 {d - 4)(d - 5f(d + 2) - 2(d - l){d - 3) 2 {d - 4)(d - 5) 3 (1 - a ) 



A 



2:-! 



-(d-5) 2 , 



A 



24 



-(d - 3)(d - 5) 2 . 



(5.11) 



Here again the coefficient of C^Ii is not gauge invariant. The on-shell matrix element (| 
is gauge-invariant; this is a strong check of our calculation. 



5.3 Renormalized decoupling coefficient 



Now we calculate the renormalized decoupling coefficient (5^). To this end, we re-express 
a s (/j,) in Zq via g 2 ,; a' s (lj) in Z'q via g^, which is, in its turn, re-expressed via g$ using 
the bare decoupling coefficient C^; combine all factors in ( |5.2[ ); and, finally, re-express the 
result via a s (/x) and expand the coefficients in e. We arrive at 



Q^) = 1 + C F T F [L 



5 r 89 

-L H 

3 36 



4-7T 



2 1 
+ C F T F l - [IIC A - 4T F (m + 2)] L 3 + - [I6(vr 2 - 6)C F - {4tt 2 - 39)C A + 20T F ] L 2 



+ 



1 



112 
~9" 



173 



-7T 



6 



C F - 76C 3 



28 



401 



-7T 



C A 



2 ( 124 
-T F 53m H 

3 V 3 



+ 8S4 



86 4 1427 

7T H 

405 27 



, n 43 4 1471 

4£ 4 7T 4 + 



1600 
243 
400 2 



-7T 



7219\ 
162 y 
3845 A 



Cs + 243^ + IST J ^ 

1 A 1B> 1685 
Tpni + - H2C 3 



81 54 " J 243 

1327 



27 



Air 



(5.12) 
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Here L is given by (|4. 15 ) and B4 by (|3.3| ). This is our main result. The two- loop part 
agrees with ||. 

At the normalization scale fj, c defined as the root of the equation 



m c (fi c ) = ii c 
we have, for the physical SU(3) colour group, 



(5.13) 



Wc) 1 + 864 I n 



+ 



1 



-2Li 4 -^ln 2 2(ln 2 2-^ 2 ) + 



427 4 815 A 175 2 877 

7T 4 C3 H 

3240 1728 s 486 108 



I if _ 1^1 
3 V ~ "86T 1 " / 



7T 



1 + 0.1030 



7T 



+ (0.7828 + 0.0124n z ) 



IT 



(5.14) 



where the number of light flavours ni = 3 includes neither b nor c. 

It is easy to re- write ( |5.12| ) via the on-shell mass m° s instead of the MS mass m c (fj,). 
Substituting the well-known one-loop relation between them, we obtain 

89^ f ajjj) 
4ir 



Cj(fi) = 1 + CpTp L 2 os - -L os + — 



36 

+ C F T F \ ~ [IIC a - 4T F (n, + 2)] L 3 os + ~ [4(4vr 2 + 3)C F - (4ir 2 - 39)C A + 20T F ] L 2 S 



88C3 
2 



112 



65 



6 



" 7r " ~ I Cp — I 76C 3 



28 . 2 -^IC 



124 

3 2>l53n, + — 



+ ( 8 ^- 405^ + 27 



1427 1600 2 9379 \ 



243 



4B 4 - «„ 



1471 



32C 3 



81 54 
1327 



400 2 

-C3 H 7T 2 + 

243 



162 

3845 \ 



486 



C A 



27 



where 



r F n; + - ( 112C 3 



L ns = 2 In 



1685 



JL 



If 



47T 



+ 



(5.15) 
(5.16) 



At \i = m° s we have, for the physical SU(3) colour group, 



sjlmc j " i + 864 



7T 



427 4 815 , 175 2 1057 

7T 4 C3 H 7T 

3240 1728 486 108 
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1 / 1327 \ 



a s (m° s ) 



TT 



3 



« 1 + 0.1030 ( ) + (0-5976 + 0.0124n;) f ) • (5.17) 

6. Conclusion 

The results on decoupling c-quark loops in the 6-quark HQET obtained in the present paper 
can be used to improve the accuracy of extracting fg from lattice HQET simulations. It 
requires the following steps (at least conceptually; in practice, some of them can be grouped 
together): 

• Matching the lattice HQET to the continuum HQET in the MS scheme at a low scale 
/i ~ l/o, where a is the lattice spacing. It can be done using the lattice perturbation 
theory, or non-perturbatively. 

• Running up to [i ~ m c in the low-energy HQET with 3 active flavours. 

• Matching the low-energy HQET to HQET with 4 active flavours. This is the subject 
of the present paper. 

• Running up to [i ~ in HQET with 4 active flavours. 

• Matching to QCD with 5 active flavours, including b. 

Currently, all steps (except the first one, which is lattice-specific) can be done at the next- 
to-next-to- leading order: the QCD/HQET matching is known at two loops ||, and 
the HQET heavy-light current anomalous dimension (which determines running in both 
HQETs) — at three loops p4|]. The present paper is a first step towards the N 3 L order 
calculation: the step 3 can now be done with the three-loop accuracy. To complete this 
program, we need the QCD/HQET matching coefficients at three loops and the anoma- 
lous dimension at four loops. The first task does not seem impossible: the QCD/HQET 



matching at three loops can be calculated using the methods of [23]. Prospects of obtaining 
the four-loop anomalous dimension are doubtful. However, the intervals of [i (from 1/a to 
m c and from m c to m&) are not really wide. Though formally the N 3 L order calculation 
requires four-loop running, its effect is likely to be small, and three-loop running should be 
sufficient. 



We have also obtained the generic formula ( A. 15| ) for the three- loop decoupling co 



efficient for the flavour- nonsinglet QCD current with n antisymmetrized 7-matrices. The 
expressions for n = 0, 1, 3, 4 can be obtained from published results; the formula for n = 2 
is new. 



The results of the present paper in a computer-readable form can be found at [25]. 

The method of calculation of three-loop on-shell HQET diagrams with massive-quark 
loops can be applied to other physical problems. We hope to consider them in a future 
publication. 
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A. Decoupling in QCD 

A.l Decoupling for the light-quark field 

In order to find the decoupling coefficient for the heavy-light current, we need the decou- 
pling coefficients for both the heavy-quark field (Sect. ^|) and the light-quark one. The 
last coefficient has been calculated in however, only the result for N c = 3 has been 
presented in the paper. Here we re-calculate this quantity. At N c = 3, our result coincides 
with Q 3 . 

Similarly to Sect. ||, the bare decoupling coefficient for the light-quark field is given by 

Z os 



K = - ztttt , Z>° s = ——— = 1, (A.2) 



where 

70s 

Jq ~ 1 - Sy(0) ' " q ~ I- S' y (0) 

and the massless-quark self-energy is S(p) = |&y(p 2 ). The renormalized decoupling coef- 
ficient is 

, Z'(a'(u),a'(u)) n 

M= 7) }/: ;,!! (° q - a.3 

Z q (a s {fj,),a{n)) 

Only diagrams with a massive quark loop contribute to Sy(0); they can be obtained 
from Fig. [5] by replacing the HQET line by a massless quark line. We made minimal 
replacements in the code calculating the HQET self-energy, so that the correctness of the 
light-quark result provides an additional check of the HQET calculation. The result is 

V ffn r T 9^ T2 (d-i)(d-2)(d-4)(d-6) 
E V (0) = C F T F ^-^I, 2d{d _ 5){d _ 7) 

-c^^& d -^t £ aw-, (a.4) 

^ ' i=l j=l,2,5 

where all non-zero coefficients are 

0"ll=Sll, CT15 = -2<T 2 5 = 815 , 0"32 = S32, &41 = «41 , CT45 = S45 , 

a 12 = 2(d - 2)(d - 3)(d - 5) 2 (4d 5 - 83d 4 + 673d 3 - 2646d 2 + 4952d - 3368) , 

cr 21 = 3(d - 5)(d - 7)(d - 10) (5d 7 - 156d 6 + 2040d 5 - 14470d 4 + 59897d 3 - 143640d 2 



3 K.G. Chetyrkin has kindly provided an unpublished three-loop formula for C?(m) f° r the SU(N C ) colour 
group; our result coincides with it. 
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+ 182340c? - 92176) , 
<7 22 = -(d-S)(d- 5) (2d 7 - 53d 6 + 588d 5 - 3557c? 4 + 12727d 3 
- 14932) - 2d(d - l)(d - 3) 2 (d - 5) 3 (1 - o ) 



26983d 2 + 31224d 



(A.5) 



(see ( |4.13D ). The on-shell renormalization constant of the light-quark field ( |A.2| ) is not 
gauge-invariant, starting from three loops; the ao-dependent term in it is the same as in 
the on-shell renormalization constant of the HQET field (4.6), see fl4.10p and ( [4.13 ). 



Now we calculate the renormalized decoupling coefficient (A. 3): 



C,(m) = 1 + C F T F [2L - 



«s(m) 



47T 



+ C F T F l -C A a(fx)L 3 



2C F 



ya( M ) + l)c A - It F 



L 



+ 



15Cf + 1 ("lT a(/i) + 232 ) Ca " ¥ 



T F ni 



L 



2387 \ 
IT) 



o(m) - 8Cs 



1187 
81 



70 

-T F (2^ + 1) 



4tt 



+ 



(A.6) 



A. 2 Decoupling for light light currents 

Here we shall consider decoupling for the flavour non-singlet QCD currents with n anti- 
symmetrized 7-matrices 



j = qoTrq , V = 7^ • • • ^ , 



(A.7) 



where r is a flavour matrix with Trr = 0, for an arbitrary n at three loops (the two- loop 
result has been obtained in (5|). 

The proper vertex F(p,p') = T + A{p,p') at p = p' = has the structure T(0, 0) = T n r, 
where V is the Dirac matrix in (|A.7|) , and T n = 1 + A n is scalar. Similarly to Sect. B the 
bare decoupling coefficient is 

£ = z° s r n , (A.8) 

because Z'° s = 1, T' n = 1. This is nothing but the on-shell matrix element of the current 
in the full theory, and it must be gauge invariant. The renormalized decoupling coefficient 
is 



, Sn J 



(A.9) 



Z n (a s (n)) 

where jo = Z n (a s (fj,))j(fi), and Z n can be reconstructed from the anomalous dimension 
j n , which has been calculated at three loops for a generic n in [26j ]. 

The diagrams for A( 3 )(0,0) can be obtained from the HQET vertex diagrams (Sect. |5.2[) 
by replacing the HQET line by a massless quark line. To find the result for an arbitrary 



- 20 - 



n, we follow the method used in [ 10 1 . If we make no assumptions about the properties of 
the matrix V, then the three- loop vertex has the structure 



A< 3 )(0, 0) = X1 T + x 27 ^V 2 r 7M2 7 m + x 3 7 M WV 4 r7/M 7^7/^1 , 



(A.10) 



because the number of 7-matrices to the left of V and to the right of it has to be even (for the 
massless quark) . Taking separately traces of 7-matrices to the left of T with Li , and of those 
to the right with R i} where U x B4 = 1 x 1, 7^7^ x 7^7^, Y 1 Y 2 Y 3 l" 4 x 7^7^ 3 7^7^, 
and solving the linear system, we can find Xi via these double traces. We may apply 
this procedure to the integrands of the diagrams, thus expressing Xi via scalar Feynman 
integrals. Now we specialize to T being the antisymmetrized product of n 7-matrices 



and obtain 



7^7^ ^7/^ =hT, h = {d-2nf 



A$ = xi + x 2 h + x 3 h 2 . 



(A.11) 



(A.12) 



We made minimal replacements in the code calculating the HQET vertex. For the 
vector current (n = 1), the Ward identity Ai = — Sy(0) gives = 1, to all orders of 
perturbation theory. Therefore, corrections in £® ( |A.8| ) are proportional to h — (d — 2) 2 . 
We obtain 



I + c h r F 4^[h-(d-2f 



+ 



9o m cO £ 



(47r) d 
1 



2 (d-2)(d-6) 
2d(d - 5)(d - 7) 



(4vr) d /2 d{d 
where all non-zero coefficients are 



1)2 2 C i I j v ij + --- 



i=l i=l,2,5 



(A.13) 



V32 = 032 
1 

Vl2 = ~ 



Vl5 = 0"l5 
V41 = 0"41 



V21 = O21 
t>45 = <745 



^25 = 0"25 



(d - 3)(d - 5) 2 (d - l)(d - 4) 2 (d - 5)(d - 9)/i 



- 4 (4d 6 - 94d 5 + 882d 4 - 4209d 3 + 10681d 2 - 13532d + 6736) 
v 22 = -(d - 3)(d - 5) [(d - l)(d - 4) 2 (d - 5) 2 /i 

+ 2d 7 - 55d 6 + 629d 5 - 3878d 4 + 13904d 3 - 28928d 2 + 32274d - 14932 



(A.14) 



The result is gauge invariant, as expected. 



Now we calculate the renormalized decoupling coefficient (A. 9): 



CM = l + -C f T f (n-l) 



o 2 1 
2(n - 3)L 2 - -(n - 15) L + — (85n - 267) 



18 



47T 



+ CpTp(n - 1)1 — (n - 3) [11CU - 4T F (n z + 2)] L 3 
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2 

+ 3 



+ 



(n - 3)(5ra 2 - 20n - 8)C F - ~(n - 3)(6n 2 - 24n - 29)C A + ^(n - 15)I> 

3 9 

16(n - 3)C 3 + ^( 17rj3 + n 2 - 275n + 117) J C* F 



16(n - 3)C 3 + — (18n 3 + 306n 2 - 2003n + 405) ) C A 
81 



4 , , 496 , 

— (163n - 477)T F m - — (n - 3)T F 
ol ol 



+ 



16(n - 3)£ 4 - -(n - 3)vr 4 + -(187n - 513)C 3 
5 3 



+ — (1727n 3 - 11681n 2 + 15946n + 12294) C F 
54 J 

- 8(n - 3)B A - -(n - 3)vr 4 + -(725n - 1887)C 3 
3 [ 5 9 

+ ^( 3087n3 - 21393n 2 + 22565n + 49881) 



C A 



+ 



4 

27 
2 

27 



32(n-3)Cs- 2^(H71n-3981) 



112(n - 3)C 3 - — (1841n - 5055) 

— / 



T F ni 



4tt 



+ 



L 



(A.15) 



The vector current does not renormalize (Z\ = 1, Z[ = 1), and we have ClO-O = 1 to 
all orders. This is to be expected; for example, for a diagonal flavour matrix r the integral 
of the vector current is a combination of the differences of the full numbers of quarks and 
antiquarks of several flavours, and these differences are some integers which are the same 
in both theories. Decoupling of the scalar current is related to that of the mass: 



Co(A0 = CmV) > 



(A.16) 



to all orders. Our result ( [A.15 ) at n = reproduces the three-loop mass decoupling (p. 

The currents j{jx) with n = 4, 3 differ from those with n = 0, 1 by insertion of the 
't Hooft-Veltman They differ from the corresponding currents with the anticommut- 

ing 7^ c by finite renormalizations: 



(A.17) 



The finite renormalization constants Zp, Z A are known with three-loop accuracy 27, 28 
Inserting 7^ c does not change the decoupling coefficient. Therefore, 



U») _ Zp(pM1 C 3 (m) _ Z' A (a' s (fi)) 

Co(m) Zp(a 8 (fj,)) ' Ci(m) Z A {a s (n)) ' 

We have checked that our result ( |A.15| ) satisfies these strong checks. 



(A.18) 
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The result for the tensor current is new: 



C 2 (/,) = i + ^t f (-2L 2 + |l-^ 



Air 



+ C F T F < 



1 

+ 9 



— [UC A - 4T F (m + 2)] L 3 + | ^28C F - ^C A - ^T F \ L 2 
-(144C 3 + 293)C F + (l44C 3 + ^p) + ^T F n, + ^T F 



L 



1 4 4 278 ^ 5639 \ 1 / _ 4 4 437 A 34135\ , 

- ( 16B 4 - 7^ + — Cs - -^r J C F + -[8B 4 - — 7r + — Cs - ttTo - ) C A 



3 V 
+ |(32 C 3 



5 ' 3 
1639 



27 



, 7>rai I 112C 3 

27 J 27 1 



3 

1373 
~27~ 
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